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1. Introduction 

Determinants appear naturally in physics when one studies systems of free fermions: 
their wave functions (Slater wave function), as well as their grand canonical partition 
function, can be expressed as determinants. These statements have exact discrete coun- 
terparts. In particular, discrete dynamics can be formulated in terms of transfer matrices 
(instead of a Hamiltonian) : these are familiar objects of statistical mechanics, which from 
a combinatorial point of view simply count the number of ways to go from a given initial 
configuration to a given final configuration. Here we are more specifically interested in 
models on two-dimensional lattices, in which the role of free fermions is played by non- 
intersecting paths; the analogous deter minantal expressions can then be derived from the 
so-called Lindstrom-Gessel-Viennot (LGV) formula [Q. In this paper we intend to show 
how these determinantal techniques can be applied to some combinatorial problems, which 
all amount to the enumeration of certain types of rhombus tilings. 

The methods we present here are fairly general and make use of various standard 
combinatorial objects. Young diagrams appear as a way to encode locations of paths 
crossing a given straight line; furthermore, the determinants involved often turn out to be 
Schur functions sy(x) where, beside the Young diagram Y, appears the set of parameters 
x = {xi, X2, ■ ■ • } which encodes the dynamics (here we are mostly concerned with simple 
counting, in which case Xi = 1). Also note that in the context of free fermions, Schur 
functions are natural building blocks for tau-functions of classical integrable hierarchies, 
see for example || for such a connection and further relations to matrix integrals. 

The paper is organized as follows. In section 2, we introduce the basic methods and 
tools which are required for our computations - the LGV formula, and the definition of 
the elementary transfer matrices from which all can be built - and revisit as an example 
the MacMahon formula. In section 3, we provide some further examples of enumeration of 
rhombus tilings of various domains. In section 4, we show how similar considerations also 
enable us to count fully packed loop (FPL) configurations with given connectivities and 
various symmetries. We conclude in section 5 with some open issues, including a discussion 
of the asymptotic enumeration of FPL configurations. 
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2. Basic ingredients: corner transfer matrices and the LGV formula 

2.1. The LGV formula 

Here and in the following we consider lattice paths whose oriented steps connect 
nearest neighboring vertices of the two-dimensional integer lattice Z 2 , and with only two 
directions allowed say along the vectors (—1,0) and (0,1). The LGV formula allows to 
express the number V({S}, {E}) of configurations of n non-intersecting lattice paths with, 
say, starting points Si, S2, ■ ■ ■ , S n and endpoints E±, E%, . . . , E n in Z 2 , as the determinant 
of the matrix whose element (i, j) is the number V(Si, Ej) of lattice paths from Si to Ej, 
namely 

V({S}, {E}) = det (V(Si, Ej)) ^. ^ N . (2.1) 

For reasons that will become obvious later on, the matrix with entries V(Si,Ej) will be 
called transfer matrix for the corresponding path counting problem. 

2.2. The matrices W and T 




Fig. 1: The transfer matrix T is expressed as a product WW 1 . 

In this note we will mainly be dealing with lattice paths with starting and endpoints 
occupying consecutive positions along two lines that intersect. The most basic situations 
are depicted in Fig. 1. The most fundamental situation is that of (ordered) starting points 
S a = {iai 1), a = 1, 2, . . . , N, and endpoints E a = (r a , r a ), a = 1, 2, . . . , N, namely with 
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paths across a corner of angle 45°. The corresponding "corner" transfer matrix W has the 
entries 

W <>P = P((i,l),(r,r))= (2.2) 

expressing that r — 1 vertical steps must be chosen among a total of £ — 1. The action of 
W is represented in the lower corner of Fig. 1. The LGV formula allows us to write 

V({(i a , l)}a=i,...,iv, {(r b , r b )} b=1 _ N ) = det (W^J^ ^ . (2.3) 

Upon reflecting the picture and exchanging the roles of starting and endpoints, we easily 
find that W l is the corner transfer matrix for the situation with starting points (r a , r a ) and 
endpoints (1, zj,). We deduce that the corner transfer matrix for the situation S a = (£ a , 1) 
and E b = (l,jb) is nothing else than T = WW 1 . We easily find that 

T hJ = V((i, 1), (1, j)) = Wi,rW j>r =( l+ ^_~ 2 ) (2-4) 

which eventually expresses that i — 1 horizontal steps must be taken among a total of 
i + j — 2. Similarly, the LGV formula gives 

^({(«a,l)}a=i,...,jv,{(l,ib)}6=i,...,jv) = det (T la Jb ) ^ ^ N . (2.5) 

Note that in both (|2.3| ) and ( |2.5| ) the result is simply the minor determinant corre- 
sponding to a specific choice of N rows and columns of the matrices W or T, taken of 
sufficiently large size. 

2. 3. Rhombus tiling of a hexagon: the MacMahon formula revisited 

As a first application of the above, let us evaluate the total number of rhombus tilings 
of a hexagon of size a x b x c, using the three elementary rhombi of size lxl. 

The standard approach consists of introducing so-called De Bruijn lines that follow 
chains of consecutive rhombi of two of the three types. In Fig. 2, we have represented 
the a De Bruijn lines connecting consecutive points along the two sides of length a of the 
hexagon. Upon slightly deforming the rhombi, as indicated in Fig. 2, these lines form a set 
of a non-intersecting lattice paths, with starting points Si = (c + £ — 1, £), i = 1, 2, . . . , a 
and endpoints Ej = (j, b+j — 1), j = 1, 2, . . . , a. Moreover, rhombus tilings of the hexagon 
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Fig. 2: A sample rhombus tiling of a hexagon of size a x b x c (left). We have 
indicated in red the De Bruijn lines following the two types of rhombi shown in the 
oval. We have deformed them so as to obtain a set of a non-intersecting lattice 
paths (right). 

are in bijection with such non-intersecting lattice path configurations. The total number 
of such configurations reads, according to the LGV formula: 

N(a, 6, c) = det (H^a)^)^^ (2.6) 

where we have introduced the "parallel" transfer matrix Hf,^ of size a x a, with entries 

Note that Eq. ( |2.6| ) expresses N(a, 6, c) as the Schur function associated to the b x a 
rectangular Young diagram and to parameters x\ = ■ ■ ■ = Xb+ C = 1 (i-e. dimension as a 
GL(b + c) representation). Similar occurrences will be discussed in more detail below. 

By line and column manipulations, this determinant may be explicitly computed to 
yield the celebrated MacMahon formula 

i=l i=l fc=l J 

We may alternatively decompose the hexagon into three parallelograms (there are 
exactly two ways of doing this, pick one), meeting at a "central point" (see Fig. 3). We now 
follow particular De Bruijn lines joining the inner edges of these parallellograms without 
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Fig. 3: Using T to enumerate rhombus tilings of a hexagon with a given number of 
winding De Bruijn loops. The first picture shows a sample rhombus tiling, and its 
decomposition into three parallelograms. In each parallelogram, the De Bruijn lines 
follow sequences of the two rhombi displayed in the ovals. In the second picture, 
we only retain the De Bruijn lines: within each parallelogram, these are clearly 
generated by the corner transfer matrix T with the appropriate shape. 

ever exiting the hexagon. These are sequences of two of the three tiles, the pairs differing 
in the three parallellograms. These lines actually form loops that wind around the central 
point. Note finally that the problem of counting configurations of lines with fixed ends 
within each of the three parallelograms uses a corner transfer matrix T of appropriate 
size. In Fig. 3, we have denoted by T(x,y) the rectangular matrix with entries ( |2.4j) , 
for % = 1,2, ... ,x and j = 1, 2, . . . , y. It is now easy to write the number N^(a, b, c) of 
configurations of exactly d non-intersecting loops winding around the central point: 

N d (a,b,c) = 1< det <d (T ar , bj T bj)Cfe T Cfe)aa ) 

1 < oi < • • - < ad < a 
1 < b t < ■ ■ ■ < b d <b 

1 < Cl < • • • < Cd < c 

l<oi < • "<ad <a 

= det(J + (iT(a, b)T(b, c)T(c, a))Ld 

aXa 

where the last line uses the decomposition formula for the determinant of the sum of two 
matrices in terms of their minors, and we have to extract the \x d coefficient to get the sum 
over all diagonal minors of size d x d. This gives a refinement of the MacMahon formula. 



(2.9) 



We recover the total number of tilings of the hexagon by summing over d, or equivalently 
by taking = 1 in the last determinant, namely 

min(a,fa,c) 

N(a,b,c)= N d (a,b,c) = det(I + T(a,b)T(b,c)T(c,a)) . (2.10) 

d=0 

Note that the symmetry under the permutations of a, 6, c is manifest here, as opposed to 
eq. (|2.6|). In the case a = b = c, denoting T(a) = T(a, a), we have 

N(a, a, a) = det(J + T(a) 3 ) = det(J + T(a)) det(J + uT(a)) det(I + u 2 T(a)) (2.11) 

where u = e 2l7r ^ 3 . Nice formulas |3|,^| happen to be known for the characteristic polynomial 
of T(a) precisely at sixth roots of unity, and allow for yet another expression for the 
MacMahon formula in a cube, via ( |2.11| ). 




Fig. 4: Decomposition of the transfer matrix Hb, c (a) for De Bruijn lines of a 
a x b x c hexagon which allows to keep track of paths passing to the left or right of 
the "center" (red dot). We give an explicit representation of the matrices W and 
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It easy to relate this picture to the former, with a set of a standard De Bruijn lines 
connecting the opposite sides of length a in the hexagon. Indeed consider the intersection 
of the standard De Bruijn lines, as well of winding lines, with the line parallel to the a sides 
that goes through the center: the location of the lines is precisely the same to the right of 
the center, whereas it is complementary to the left. Therefore Nd(a,b,c) also counts the 
tiling configurations with exactly d of these lines passing to the right of the "center" of the 
hexagon. This is also a consequence of the following relation between transfer matrices, 
pictorially proved in Fig. 4: 

H b , c (a) = W(b, ayP b T(b, c)P c W(c, a) + Z h {a)Z c {af (2.12) 

in which the [if) entry of the first term counts the number of paths passing to the right 
of the center connecting the i-th entry point to the j-th exit one, while the second term 
yields the contribution of paths passing to the left of the center. In ( |2.12| ), as usual the 



argument corresponds to the size of the matrix ((a) for a square a x a matrix, (a, b) for a 
rectangular a x b matrix), and the matrices have the following entries, see Fig. 4: 

(Hb, c )ij = + { ._ A , Wij = ( _ A , [P b )ij = Sj t b+i-u {Zb)ij = C_j 

(2.13) 

Note that the matrix P induces a reflection, needed whenever passing from a corner to 
another one pointing in the opposite direction. We finally identify 

det(I + fiT(a,b)T(b,c)T(c,a)) = det (Z 6 (a)Z c (a)* + pV{b, afP b T{b, c)P c V{c, a)) (2.14) 

a X a a X a 

by noticing that, after factoring out ZZ l whose determinant is 1, the two corresponding 
matrices are conjugate of one-another. 



3. Rhombus tilings of various domains 

3.1. Lozenge with glued sides 

T may be used to compute the numbers N^{a) of tilings of a lozenge of side a glued 
along two of its consecutive edges (see Fig. 5), and on which exactly d De Bruijn loops 
wind around the tip of the cone. We simply have 

N^(a) = det(I + ^T(a))\ fld . (3.1) 
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Fig. 5: Computing the numbers of tilings of a lozenge of side a glued along two of 
its consecutive edges, and with exactly d loops winding around its conic singularity. 



Disregarding the d dependence, we find the numbers 

N L (°) = ^I + T(a))=Yl \ * > . (3.2) 

n=0 V ' V ' 

These read for a = 1, 2, 3, . . . 

2, 5, 20, 132, 1452, 26741, 826540, 42939620, . . . (3.3) 

This is by construction the number of cyclically symmetric plane partitions in a cube 
ax ax a, obtained by decomposing a regular hexagon of size ax ax a into three identically 
tiled parallelograms, with suitable identifications of boundaries (see section 4.3 for another 
interpretation of these numbers). 

3.2. Half-hexagon with glued sides 

The corner transfer matrix W may be used to compute the numbers N^ H (a) of tilings 
of a half- hexagon of side a glued along two of its cut edges (see Fig. 6), and on which exactly 
d De Bruijn loops wind around the tip of the cone. We have 

iVf H (a) = det (i + n W(a)W(aYW(a))\, d . (3.4) 

aXa 

Disregarding the d dependence, we find the numbers 



N HH (a) = det (/ + WW l W) 



nn T-rn-2 (4j+3)!(j + l)!j! A (4n-l)!n!(n-l)! ., _ „ 

Z Ull=0 727+W7727+IW I (2n)\((2n-lV.Y 2 ^ a — £11 



Lj=0 (2j+2)!((2j + l)!)* J (2»)!((2n-l)!) a 

on+1 ( n n-l (4j+3)!(j + l)!j! \ 2 -r _ „ , -, 

^ ^llj=0 (2j+2)!((2j + l)!)^ J 

(3.5) 
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Fig. 6: How to use the corner transfer matrices W and W t to generate the rhombus 
tilings of a half- hexagon with glued sides. The gluing is indicated by the arrows. 
The half-hexagon is decomposed into three triangles, within which the De Bruijn 
lines are generated by W or W f . 

These numbers read for a = 1, 2, 3, . . . 



Remarkably, these turn out to match exactly the total dimension of the homology of free 
2-step nilpotent Lie algebras of rank a || (cf entry A078973 of the on-line encyclopedia of 
integer sequences ||). 

Let us give an explanation for this apparent coincidence, and present a determinant 
formula for the Poincare polynomial of 2-step nilpotent Lie algebras of rank a. As in the 
case of the full hexagon, one can introduce the standard De Bruijn lines; here there are a 
lines which enter from the middle side of length a and come out from the side of length 
2a. Let us denote by 2a > ri\ > ■ ■ ■ > n a > the locations of their endpoints. Then the 
LGV formula tells us that the number of such De Bruijn lines at fixed endpoints nj is: 



Now the irreducible representations of GL(a) are indexed by Young diagrams Y = 
{Ai, A2, • • • , A a } with Xi boxes in the i-th row, and Aj > Aj+i > for i = 1, 2, . . . , a — 
1. The corresponding character, evaluated on the class of matrices with eigenvalues 



2, 6, 36, 420, 9800, 4527650, 41835024, 7691667984 . . . 



(3.6) 




(3.7) 
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x = {xi,X2, ■ ■ ■ ,x a } is the Schur function, expressed via the Jacobi-Trudi formula as 
s Y (x) = det (hx i+ j-i(x)) 1<iJ<a = det (^+i-<( x )) i< < j< where ^mfc) are the com- 
plete symmetric functions, generated by ^2 rn > h rn (x)t Tn = ]^[" =1 1/(1 — taj), similarly 
S m >o ^m( x )^ m = rii=i(l + ^i); an d the denote the A's of the transposed Young dia- 
gram Y T (reflected w.r.t. the first diagonal). The dimension dimy of the representation 
Y correponds to the character of the identity class, with x\ = x% = ■ ■ ■ = 1 (also denoted 
by x = 1). Explicitly, 



dimy = sy(l) = det 



a + Xi+j — i — 1 
K+j-i 



det 



l<i,j'<a 



A !+J-Vi<i, ? 



(3.8) 



We recognize in Eq. (|3.7|) this dimension, provided we identify rn = \[ + a + 1 — i. Finally, 
the gluing of the two half-sides means in the summation over the endpoints rii, we must 
only include those such that the corresponding diagram Y is invariant by transposition. 
We thus find: 

N HH {a) = dim Y ( 3 - 9 ) 

Y=Y T 

which is precisely the quantity calculated in M. 




m. 



m 



Pi 



P 2 



Pc 



(b) 



m. 



Fig. 7: (a) Half-hexagon tilings with d De Bruijn lines having fixed endpoints, 
with positions mi, m2, ■ . . , md and pi,p2, ■ ■ ■ ,Pd counted from the center, (b) The 
corresponding Young diagrams, with characteristics (m;p) in Frobenius notation. 
Boxes are counted from and include the diagonal, from left to right (m's) and top 
to bottom (p's). 
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Note that the correspondence between the two approaches can be made for each term 
of the sum: indeed the endpoints of the standard lines and of the winding lines are clearly 
in one-to-one correspondence. If we record the positions of endpoints of the winding lines 
on the (cut) boundary, counted from the center, namely (m;p) = {mi, . . . , m^pi, ■ ■ ■ ,Pd}, 
1 < mi < 7B2 < • • • < ma < a from center to top and 1 < p\ < P2 < ■ ■ ■ < Pd < a from 
center to bottom in the case of d winding lines, then the corresponding Young diagram Y 
is given by Fig. 7. Therefore, the total number of tiling configurations N^ p ^(a) of the 
half-hexagon with fixed endpoints (m;p) of De Bruijn lines is also the dimension of Y: 

^S)(<0 = det ((^^Wp,) = dimY ■ ( 3 - 10 ) 

The Poincare polynomial P a (u) generating the dimensions of the homology spaces of 
a 2-step nilpotent Lie algebra of rank a has been shown to read 

P a (u) = dim ru ]Yl (3.11) 

Y=Y T 

where the sum extends over the self-transposed Young diagrams, namely those for which 
the two characteristics (m;p) are equal, and where \Y\ denotes the total number of boxes 
in Y. Note that |y(ra;ra)| = Xli<i<d(2 TO i — !)• Introducing the diagonal matrix 9 with 
entries 

6 ij = S ij u 2i - 1 , i,i = l,2,...,a (3.12) 
and using the relation ( p.lU| ), we get 

P a (u) = det(7 + eWW l W) . (3.13) 

a X a 

The first few such polynomials read 

Ptiu) = l + u 

P 2 (u) = l + 2u + 2u 3 + u 4 

P 3 (u) = 1 + 3u + 8u 3 + 6w 4 + 6m 5 + 8u 6 + 3u 8 + u 9 

P 4 (m) = 1 + Au + 20u 3 + 20u 4 + 36u 5 + 64u 6 + 20u 7 + 90u 8 + 20u 9 + 64u 10 + 36U 11 

+ 20u 12 + 20u 13 + 4u 15 + u 16 
p 5 ( u ) = 1 + 5u + 40u 3 + 50w 4 + 126m 5 + 280m 6 + 160m 7 + 765m 8 + 245m 9 + 1248m 10 + 720m 11 

+ 1260m 12 + 1260m 13 + 720m 14 + 1248m 15 + 245m 16 + 765m 17 + 160m 18 + 280m 19 + 126m : 

+ 50m 21 + 40m 22 + 5m 24 + m 25 . 

(3.14) 

At m = 1, we recover the formula fl3.5| ), which indeed counts the total dimension of the 
homology of a 2-step nilpotent Lie algebra of rank a. 
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b 




b+m 

Fig. 8: The transfer matrix T m (a,b) is used to generate tilings of a hexagon with 
a central triangular hole. 

3.3. Hexagon with a triangular central hole 

A variant of T may be used to generate the numbers Nd(a, c + m,b,a + m,c,b + m) 
of tilings of a hexagon with an equilateral central triangle of side m removed (see Fig. 8), 
and with d loops winding around the hole. More precisely, let T m (a, b) denote the transfer 
matrix of rectangular size a x b with entries 

, m % ffTi + i + j — 2\ 
with i = 1, 2, . . . , a, j = 1, 2, . . . , b. The configurations are then generated by 

min(a,b,c) 

det (I + iiT m {a, 6)T m (6, c)T m (c, a)) = > /j, d Nd(a,c + m,b, a + m, c, 6 + m) . (3.16) 

aXa ^— ' 

d=0 

See for explicit expressions of some of these determinants. 

3.4- Hexagon with chopped- off corners 

We consider the tiling of a hexagon axbxc with three corners chopped off as indicated 
in Fig. 9, along line cuts at distances m : p : q from the central point, with m > max(a, b), 
p > max(6, c), and q > max(a, c). Let N^ rn,p,q \a, 6, c) denote the number of tilings of this 
domain with d winding lines. These are enumerated using a restricted version of T that 
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a 




Ab,c) 



Fig. 9: A hexagon of shape a x b x c with three chopped off corners at distances 
m,p,q from the origin (central point). We have represented the relevant corner 
transfer matrix. 

incorporates a "ceiling" at a given height not to be crossed. By the reflection principle, 
the relevant truncation reads 

'i+j-2\ _ fi + j -2 N 
i — 1 / \ m 



rp{m) __ 
ij 



(3.17) 



We denote by T( m )(a, b) the corresponding rectangular matrix a x b. With this definition, 
we have 

det (i + iiT (m) (a,b)T( p \b,c)T(i\c,a)) L d = N { d m ' p ' q \a,b,c) . (3.18) 

aXa \ / 

Note that a half-corner transfer matrix can still be introduced, namely a matrix 
with entries 

A: _ 1 \ / i _ 1 \ 

(3.19) 



W {m) 



% — 1 \ / i- 1 
v i - V v m + 1 - h 

and such that T^ m \a,b) = W^ m \a, c)W^ m \b, cf for any c > m/2. 

A particular case of (|3.18| ) corresponds to tilings of a notched equilateral triangle, 
with a = b = c and m = p = q = a. 

We may easily treat the case of rectangular chopping as well, by suitably modifying 
T to incorporate the corresponding "broken ceiling" restriction indicated in Fig. 10. Let 
Nd(ai, 0,2] bi, 62; ci, C2) denote the corresponding number of tilings with d winding loops. 
We have the transfer matrix 

T{a u b 1 ) 2o ibl (ai,6 2 ) 
T ai ,o(a2, bx) T(a 2 , h)T(bi t a x )T{a x , b 2 ) 



T(a 1 ,a 2 ;b 1 ,b 2 ) 



(3.20) 
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Fig. 10: Tilings of a hexagon with removed parallellograms. The relevant transfer 
matrix T has a block decomposition according to the arrows. 

where T m , p (a, % = ( m+ ^-i" 2 )' 1 < * < a > 1 < 3 < b, and finally 

A^ d (ai,a 2 ;6i,6 2 ;ci,C2) = 

det (I + iiT(a 1 ,a<2;b 1 ,b2)T(b 1 ,b2;c 1 ,c< 2 )T(c 1 ,C2;a 1 ,a 2 )) | M d . 
(01+02) x (01+02) 



4. Fully-Packed Loops and more tiling problems 

^.i. Fully-Packed Loops, Alternating Sign Matrices and Plane Partitions 

The fully-packed loop (FPL) model plays a central role in the Razumov-Stroganov 
conjecture, which attracted a lot of attention recently. In short, the FPL model is a 
statistical model on a square grid of size a x a of the square lattice, in which edges may 
be occupied or not, and with the constraint that exactly two incident edges are occupied 
at each vertex. Moreover, one imposes alternating boundary conditions along the border 
of the grid, that every second edge at the exterior of and perpendicular to the boundary is 
occupied. These occupied external edges are labeled 1,2,..., 2a. The FPL configurations 
are in bijection with alternating sign matrices (ASM) of size a x a, namely matrices with 
entries ±1,0 only, such that 1 and —1 alternate and sum up to 1 along each row and 
column. 

The Razumov-Stroganov conjecture involves refined FPL numbers, according to the 
connectivity of external edges. Indeed, from the definition of the model, the occupied 
external edges are connected by pairs via chains of consecutive occupied edges, forming 
lines (rather than loops) on the square grid, while closed loops may occupy some of the 
inner edges of the grid. To summarize these connectivities, one generally uses the language 
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of link patterns, namely planar permutations n G S^a with only 2-cycles indicating the 
connected edges. The Razumov-Stroganov conjecture relates the numbers of FPL with 
fixed connectivities to the groundstate vector of the 0(1) loop model 

Rhombus tiling configurations of a hexagon are also called plane partitions (PP), as 
they may be interpreted as the view in perspective from the (1,1,1) direction of a piling- up 
of unit cubes in the positive octant of the 3D lattice Z 3 , with the constraint that gravity 
has the direction —(1,1,1) and that only stable configurations are retained. The total 
number of FPL or ASM of size a x a matches that of so-called totally symmetric self- 
complementary plane partitions (TSSCPP) of size 2a x 2a x 2a, namely plane partitions 
being maximally symmetric, i.e. under rotations of 2n/3 and reflections w.r.t. medians, 
as well as identical to their complement. This result is one of the keystones of modern 
combinatorics and is beautifully described in the book [[| . From this we simply retain that 
there should exist a natural bijection between the two sets of objects, still to be found to 
this day. 

4-2. FPL with fixed sets of nested lines 

In parallel to the FPL-ASM-TSSCPP relation, there exist bijections between special 
FPL configurations and rhombus tilings of special domains. The simplest of these concerns 
the FPL of size nxn with three sets of nested lines, in numbers say a, 6, c with a + b+c = n, 
namely with a link pattern that connects 2a external edges say 1, 2, . . . , 2a by symmetric 
pairs (1, 2a)(2, 2a — 1), . . . , (a, a + 1), then the 2b next analogously and the 2c remaining 
analogously, thus forming three "bundles" of respectively a, 6, c connecting lines. These 



were shown to be in bijection with the rhombus tilings of a hexagon of size a x b x c (T(J 
and henceforth are enumerated by the MacMahon formula 



This result was improved in so as to include the case of four sets of nested lines, in 
numbers say a, 6, c, d. The total number A^(a, 6, c, d) of such FPL matches that of rhombus 
tilings of the domain of Fig. 11, which have exactly d De Bruijn loops crossing the gluing 
line of length c + d. This decomposition allows one to write the number as 

N(a,b,c,d) = det (I + T(b + c,a + d)T(a + d,c + d + e)U(c + d,b + c))L d (4.1) 

(&+c)x(6+c) 

where the transfer matrix U(p, q) has the block form (see the right sketch of Fig. 11) 



U(p,q) = (T(p,q-p) fiX q . p {p)) 
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(4.2) 



c+d 




c+d 




b+c 



U(c+d,b+c) 



Fig. 11: The domain whose rhombus tilings match the N(a, b, c, d) FPL with four 
sets of nested lines a, b, c, d. It is a heptagonal domain with two sides glued, and 
with the indicated size. Moreover, it must have exactly d De Bruijn loops (thick 
lines) crossing the gluing line. The transfer matrices are obtained by decomposing 
the domain into three parallelograms as indicated. 




Fig. 12: FPL configurations with four sets of nested arches separated into two 
subsets by a fifth one. 

where {X m ) itj = C^ 1 -/"') • 

Finally, it is possible to extend this bijection to FPL with link patterns with five sets 
of nested lines as depicted on fig. 12. These FPL are easily related to the rhombus tilings 
of the domain displayed in Fig. 13, which have exactly d De Bruijn loops going across the 
gluing line. The only difference with the case of four nested sets of lines is that we have 
introduced a branch cut of length e in the glued domain. The numbers N(a,b\e\c,d) of 
desired FPL reads 

N(a,b\e\c,d) = det (l+T(b+c, a+d)T(a+d, c+d+e)U (e) (c+d+e, b+c)) L* (4.3) 

(b+c) X (b+c) 

where the transfer matrix has the block form 

(p, q) = I T(p, q - p + e 



q—p+2e 



(p-e) 



(4.4) 
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with a zero block of size ex (p — e). Note that all of the currently known numbers of FPL 
configurations with prescribed connectivities (for which a proof exists), as in [10JTT,12| , 
are special cases of our pattern (a, b\e\c, d). 




Fig. 13: The domain whose rhombus tilings match the N(a, b\e\c, d) FPL with a, b 
and c, d nested lines separated by e. The nonagon is glued along its edges of length 
c + d as indicated by blue arrows. The red segments of length e are forbidden, and 
form a branch cut in the glued domain. After decomposing the domain into three 
sub-domains, we have indicated the transfer matrices needed for the enumeration. 
As before, we must have a total of d De Bruijn loops going through the glued cut 
of length c + d. 

4-3. Other symmetry classes of FPL/ASM 

The Razumov-Stroganov has been generalized to symmetry classes of FPL: the num- 
bers of FPL with certain symmetries and with given connectivity patterns turn out to be 
related to the ground state vectors of the same 0(1) loop model, but with different bound- 
ary conditions |l3|Jl^| . The counting of FPL with connectivity (a, b\e\c, d) described above 
can be naturally restricted to symmetry classes by imposing the corresponding symmetry 
on the rhombus tiling, or equivalently by dividing out the domain by the symmetry. Here 
we consider only the symmetry classes for which a Razumov-Stroganov type conjecture is 
known. 

The first class is the so-called half-turn symmetric FPL (HTSFPL), i.e. configurations 
which are invariant by rotation of it. Naturally, the connectivity pattern itself must be 



17 



e 

J 

a+b 



Fig. 14: Domain of rhombus tilings for (a,b\e\a,b) HTSFPL. 

half-turn symmetric, so that we must choose it to be of the form (a, 6|e|6, a). The domain 
is of the form of Fig. 14, and we find that 

N HT (a,b\e\a,b) = det (I + [i T e (a + b)) (4.5) 

where T e (a + b) = T(a + b)X e (a + b) has its entries given by Eq. ( |3.15| ) with 1 < i, j < a + b. 

Note in particular that for e = we recover the sequence A045912 of || that was 
already mentioned in section 3.1 (Eq. (|3J])); and by summing at fixed n = a + b, the total 
number of dimers det(7 + T{n)) (Eq. (|3~2|); sequence A006366 of @) which is the number 
of cyclically symmetric plane partitions in the n-cube. This number also coincides with the 
ratio of numbers of 2n x 2n half-turn symmetric ASM (HTSASM) and of n x n ASM. One 
now finds an indirect proof of the latter fact, via the Razumov-Stroganov conjecture for 
HTSFPL and FPL. Indeed, there is a projection that sends the ground state eigenvector 
counting FPL onto the one for HTSFPL. All connectivity patterns that contribute to 
N HT (a : b\0\a, b) are projectedB onto the trivial pattern of n arches for FPLs of size n; 
there is only one corresponding FPL, and therefore the sum of such contributions must 
be related to the relative normalization of the ground state eigenvectors, which can be 
obtained by taking the sum of all components, i.e. which is equal to the ratio of numbers 
of 2n x 2n HTSASM and ofnxn ASM. 

1 Recall that a HT-symmetric pattern of size 2ra is projected onto a pattern of size n by cutting 
the disk into 2 equal sized pieces in any way such that it does not cut an arch, taking one half-disk 
and gluing it back together into a disk. 



a+b+e 
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a+b+e 




a+(e-1)/2 




(e-1)/2 



(i) 



(n) 



Fig. 15: Domains of rhombus tilings for the (i) (a,b\e\b, a) VSFPL and (ii) 
(a,l,a\e\a,l,a) HVSFPL configurations. In each case there are exactly a non- 
intersecting paths (only one is depicted). 

Next we consider the vertically symmetric FPL (VSFPL), which are invariant by 
reflection with respect to the vertical axis, and only exist for odd sizes. The connectivity 
is now (a, b\e\b, a), e odd, and the domain is of the form of Fig. 15 (i). There are exactly a 
non-intersecting paths, and it is convenient to redefine them as coming out of the side of 
length a and propagating to the opposite side as depicted. As in section (3.4), we apply 
the reflection principle to compute the number of paths in the presence of a "wall", and 
then apply the LGV formula: 



N v (a, b\e\b, a) 



det 

i,j = l.. 



2b + e + j-l 

b-j + i 



2b + e + j-l 
b-j-i+1 



(4.6) 



Finally, the vertically and horizontally symmetric FPL (HVSFPL) are invariant under 
reflections with respect to both axes. Here the sequence of arches must be slightly modified 
in order to accomodate the additional symmetry; with the insertion of two single arches, the 
connectivity pattern becomes (a, 1, a|e|a, 1, a) and the domain is of the form of Fig. 15 (ii). 
We finally find 



N 



HV, 



a, 1, alela, 1, a) 



det 

i,j=l...a 



2a + (e-l)/2-j + 1 
a - 2j + i + 1 



2a+(e-l)/2-j + l 
a - 2j - i + 2 



(4.7) 
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4-4- Plane partitions: q-decoration 



When we view the rhombus tilings of a hexagon as plane partitions, we may introduce 
yet another catalytic variable, namely a weight q per unit cube in the PP. Many of the 
above results may thus be "(/-decorated" whenever there exists a PP interpretation. This 
is for instance the case for the total number of PP in a box of size a x b x c, which gives 
a g-MacMahon formula. To express it, we introduce (/-deformed corner transfer matrices 
w, t that keep track of the numbers of boxes in the PP picture. This leads to 



w 



g6 T 2 T 



% - 1 

3-1 



(4- 



where we have used (/-binomials 



;i- g n-i)(i_ 9 w-2)...(i_g»)/((i _,)(!_ 



q 2 ) . . . (1 — q a b ))), and where the prefactors are ad-hoc to take care of boundaries. Simi- 
larly, we have t = ww 1 reading 



t 



1 , i+j-2 



'J 



i+j-2 
i-l 



(4.9) 



J <7 



With this definition, we get the (/-deformation of the MacMahon formula (|2.6|) 



a b 



N(a, 6, c; q) = det(J + t(a, b)t(b, c)t(c, a)) = JJ JJ JJ - 



i=lj=l fc=l 



_ q i+j+k-2 



(4.10) 



for the generating function N(a, 6, c; g) for PP in a box a x b x c with a weight g per unit 
cube. Combining this with eq.(|2.9|), we also get the generating function N(a, 6, c; q, //) for 
PP in a box a x b x c with a weight g per unit cube, and with a weight \i per De Bruijn 
loop in the tiling picture 



N(a, 6, c; g, //) = det(7 + /j,t(a, b)t(b, c)t(c, a)) . 



(4.11) 



With the matrix t, we have also access to the generating function N cs (a; g, //) for the 
numbers of cyclically symmetric plane partitions in a box ax ax a, weighted by q per unit 
cube, and \i per winding De Bruijn loop, namely 



N (a; q, (j) = det(J + n t(a)) 
20 



(4.12) 



1 




q 5 \X q 7 \X q 8 JO. 2 



Fig. 16: Vertically symmetric plane partitions in a box a x a x a corresponding to 
the tiling of a glued half-hexagon, for a = 2. For each partition, we have indicated 
the contribution to the generating function N HH (2; q, /x) = 1 + g^u + q 3 /j, + g 5 /x + 

1 8 2 

q fi + q n . 

provided q is replaced by q 3 in the definition of t (as the total number of boxes is three 
times that in any of the three identical copies of the parallelogram that form the total 
hexagon) . 

Finally, we get the generating function N HH (a; g, fi) for the numbers of vertically 
symmetric partitions corresponding to the tiling of the glued half-hexagons of section 3.2, 
with an additional weight q per unit cube and \i per winding De Bruijn loop 

N HH (a; q, y) = det (J + yw{a)w t {a)w{a)) (4.13) 

provided q is replaced by q 2 in the definition of w (as the plane partition is obtained by 
completing the (opened) half-hexagon with its reflection. These partitions are represented 
in Fig. 16 for a = 2. We may also write a (/-deformed version of the Poincare polynomial 
( pMl| ) for 2-step nilpotent Lie algebras of rank a: 

PJu; q) = det (I + —Ow(a)wHa)w{a)) (4.14) 

with the matrix 9 as in ( |3.12| ) . The latter corresponds presumably to a sum over Schur func- 
tions of the form ^2 y =y t s Y(x q )u^ Y ^ , with the specialization x = x q = {1, q, q 2 , . . . , q" 1-1 }. 
The coefficients of P a (u;q) as a function of u are g-deformed dimensions of homology 
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spaces, and await a good interpretation in the Lie algebraic context. For illustration, we 
have 

P s (u;q) = l + u(l + q + q 2 ) + u 3 q 2 (l + q) 2 (l + q 2 ) + u A q 2 (1 + q 2 )(l + q + q 2 ) 

+ u 5 q 5 (l + q 2 )(l + q + q 2 ) + u 6 q 6 (l + q) 2 (l + q 2 ) + u 8 q 10 (l + q + q 2 ) + u 9 q 12 . 

(4.15) 



5. Conclusion 

In this paper, we have provided expressions for the numbers of rhombus tilings of 
certains regions of the plane with some gluing conditions. Sometimes these expressions were 
directly determinants, sometimes they were coefficients in the expansion of a determinants 
in powers of a variable, in a way similar to grand canonical vs canonical partition functions. 
We used transfer matrix techniques, in analogy with Baxter's corner transfer matrix. We 
have applied these ideas to the computation of the numbers of FPL for a fairly general 
class of connectivities, with and without additional symmetries. Various issues remain 
open. 

5.1. FPL and rhombus tilings 

In view of sect. 4.2 and 4.3, one may wonder whether there exists a deeper connection 
between FPL configurations with prescribed connectivities and rhombus tilings of possibly 
cut and/or glued domains of the plane. 

The answer to this question must be subtle as was already observed in the case of 
four sets of nested lines [11|: Indeed in that case, only a few rotated versions (in the sense 



of Wieland |pL5| ) of the FPL counting problem allowed for a straightforward bijection with 
rhombus tilings of the domains of figure 11. In the general case, we must first optimize 
the Wieland rotation in order to attack the problem. 

Another relation to rhombus tiling, though very hypothetical, would first involve 
finding a bijection between FPL/ ASM and totally symmetric self-complementary plane 
partitions (TSSCPP), themselves reducible to the rhombus tilings of triangles with free 
boundary conditions along one side. The number of TSSCPP in a box of size 2a x 2a x 2a 
was shown to match that of a x a FPL/ ASM, but no natural bijection is available yet 
between the two sets 0. If we knew such a bijection, we could characterize among TSSCPP 
the rhombus tilings corresponding to FPL configurations with fixed connectivities, thus 
provide an answer to the above question. This answer, however, would not coincide with 
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that of sect. 4.2 in the case of four sets of lines separated by a fifth one, as the domain 
considered do not match in any simple way. 

Answering this question remains an interesting challenge. 



5.2. Asymptotic enumeration 

A by-product of our work has been to provide us with new data for FPL configuration 
numbers, either by explicit determinant formulae, or just numerical. It is tempting to 
try to extract the asymptotic behavior for large configurations of a given type. This is 



in particular the case for the configurations denoted ( in []IB|, i.e. made of p sets of r 



arches, for which we conjecture that p large, r fixed, or vice versa (and n = pr) 

log # (( )») « kt 2 {{p - l) 2 - 1) + O(logn) (5.1) 
with the same k — \ log || that appears in the asymptotics of A n , 

A n ^ n - 23/36 e Kn2 . (5.2) 



p\ 


r 1 


2 


3 


4 


2 


1 


1 


1 


1 


3 


2 


20 


980 


232848 


4 


7 


3504 


118565449 


266866085641550 


5 


42 


5100260 


1637273349805800 




6 


429 


60908609580 






7 


7436 


5939300380261111 






8 


218348 


4717858636573174999768 






(3 r 2 


~ 2? 


2 ' 





a. 





3/2 
'4 (*) 
15/2? 



Table 1. Numbers of FPL configurations ()? with p sets of r arches 



(*) see the text for justification. In each column, as p — » oo, there is good 
evidence that the numbers behave as (||) , with numbers /3 r as shown; 

a r 2 

likewise, along each row, r — >• oo, behavior as (||) p , with a p as shown. 
The expression in ( |5.1| ) is consistent with these data. 
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The conjecture ( |5.1|) agrees with some particular cases. For r = 1, all p = n, we have 
configurations with n simple arches, whose number is A n _i, as conjectured in ||, and 
whose asymptotics is thus given by ( |5.2| ) (with n \— ► n — 1). For p = 3, r = n/3, we have 
the MacMahon formula, the asymptotics of which is easy to calculate. Our conjecture is 
also supported by the observation that the number of 4-arch configurations ( )f, r = n/4, 
is given by the middle term in the expansion of det(x + T 2 ), with T the Pascal matrix 
( |2.4| ), and is in fact the dominant coefficient in that polynomial in x. Then Mitra-Nienhuis 
||17|| conjecture on det(zl + T) ~ (A HT (L = 2r)) 2 together with the known asymptotic 

L 2 /4 

behavior of half-turn symmetric ASM, A UT (L) ~ (||) , gives log#( ~ i8r 2 log||. 
Last but not least, our conjecture (|5.1|) is well supported by numerical data, as shown in 
Table 1. 
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